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POSTSWIRL  PROPULSORS  -  A  DESIGN  METHOD  AND  AN  APPLICATION 


by 

Benjamin  Y.-H.  Chen 

Carderock  Division,  Naval  Surface  Warfare  Center,  Bethesda,  MD  20084-5000,  U.S.A. 


ABSTRACT 

A  design  method  for  postswirl  (PS)  propulsors  is  ad¬ 
dressed.  The  principles  of  momentum,  mass,  and  circula¬ 
tion  conservation  are  satisfied  and  the  effects  of  the  hub 
boundaries  have  been  taken  into  account  in  the  design 
and  analysis  methods.  A  PS  propulsor  was  designed  for 
close- to- uniform  flow  at  the  operating  point  for  a  sur¬ 
face  ship.  Self-propulsion  experimental  results  show  that 
the  performance  predictions  agree  well  with  the  measure¬ 
ments. 

INTRODUCTION 

There  has  been  a  renewed  interest  in  finding  more  ef¬ 
ficient  and  quieter  propulsors  for  new  naval  surface  ships. 
A  promising  candidate  is  a  postswirl  (PS)  propulsor  which 
consists  of  a  stator  installed  behind  a  rotor.  One  of  the 
earliest  work  for  PS  propulsors  was  developed  by  Wagner 
[1]  PS  propulsors  exhibit  many  advantages  in  terms  of 
powering  and  cavitation  over  single  rotation  (SR)  pro¬ 
pellers. 

The  improvements  in  the  propulsive  efficiency  for  the 
PS  propulsor  are  due  to  the  following  reasons: 

1 .  Reduced  rotation  losses  in  the  propeller  slipstream. 

2.  Reduced  axial  kinetic  energy  losses  in  the  propeller 
slipstream. 

The  first  improvement  gain  results  from  the  fact  that 
the  tangential  velocity  generated  by  the  stator  is  oppo¬ 
site  to  that  generated  by  the  rotor  on  the  stator  plane. 
The  second  source  comes  about  because  the  stator  causes 
the  maximum  rotor  blade  loading  to  shift  inboard.  This 
only  occurs  when  the  propulsor  operates  in  a  non-uniform 
flow. 

The  improvements  in  the  cavitation  inception  speed 
for  the  PS  propulsor  result  for  the  following  reasons. 

1  Reduced  blade  loading  for  blade  surface  cavitation. 

2  Reduced  tip  circulation  for  tip  vortex  cavitation. 


The  first  item  is  due  to  the  stator  blades  taking  some 
of  the  loading  from  the  rotor  blades.  The  second  item 
is  owing  to  a  shift  of  the  maximum  rotor  blade  loading 
inboard. 

The  purpose  of  the  present  study  was  to  develop  a 
design  method  for  the  PS  propulsor  based  on  rational 
hydromechanics  using  the  principles  of  momentum,  mass, 
and  circulation  conservation.  The  method  accommodates 
non-zero  loading  at  the  blade  roots  in  the  design  and 
analysis  process  and  accounts  for  the  effects  of  the  hub 
boundaries. 

A  PS  propulsor  was  designed  for  close-to-uniform  flow 
at  the  operating  point  for  a  surface  ship.  Self-propulsion 
measurements  for  the  PS  propulsor  design  are  given. 

A  DESIGN  METHOD  FOR  POSTSWIRL 
PROPULSORS 

DESIGN  PRINCIPLE 

Momentum  conservation  requires  that  the  PS  propul¬ 
sor  generate  a  net  force  to  overcome  two  types  of  drag; 
the  bare  body  drag  and  the  drag  due  to  propulsor- hull 
interaction.  There  are  two  types  of  propulsor-hull  in¬ 
teractions;  thrust  deduction  and  wake  fraction,  but  only 
thrust  deduction  affects  momentum  conservation.  Thrust 
deduction,  an  additional  drag,  results  from  the  propul¬ 
sor  accelerating  the  local  flow  field  about  the  ship  and 
causing  a  reduction  of  the  local  hull  surface  pressure.  In 
addition,  the  wall  shear  stress  increases  due  to  a  higher 
velocity  which  in  turn  increases  the  frictional  resistance. 
The  thrust  deduction  fraction  for  propulsors  is  defined 
by 


where  Rj  is  the  bare- hull  resistance  and  T  is  the  total 
thrust  for  all  propulsors. 


1 


I  sing  Lagally’s  theorem,  Weinblum  [2]  developed  a 
procedure  to  calculate  thrust  deduction.  Cox  and  Hansen 
3,  developed  a  procedure  of  thrust  deduction  for  a  sin¬ 
gle  rotation(SR)  propeller  in  terms  of  the  potential  flow 
about  the  hull  and  appendages  which  are  represented  by 
surface  singularity  distributions.  However,  the  thrust  de¬ 
duction  factor  for  PS  propulsors  is  still  a  research  topic 
because  of  the  complex  interactions  between  the  rotor 
and  the  stator.  Since  there  is  no  thrust  deduction  mea¬ 
surement  for  the  PS  propulsor  in  the  current  study,  the 
thrust  deduction  factor  for  SR  propellers  is  employed  in 
the  present  design  method.  When  the  thrust  deduction 
value  for  a  PS  propulsor  is  measured  on  the  proper  hull 
form,  it  would  be  appropriate  to  use  the  thrust  deduction 
for  the  PS  propulsor. 

Mass  conservation  determines  the  circulation  distribu¬ 
tion  of  the  stator  once  the  circulation  distribution  of  the 
rotor  is  specified.  Based  on  mass  conservation,  the  pre¬ 
liminary  diameter  of  the  stator,  D ,,  can  be  determined 
from  the  rotor  diameter,  Dr,  through  the  following  for¬ 
mula. 


D,  =  /  x  Dr 

where 

f  —  IXm,*  (  ^Tr  ~b  Urrr  ^  rr  t)  X  Rr  d  X  Rr  I 

+  «»«  +  VzJXtodXn.  ^ 

where  VZT  and  Vr,  are  the  axial  inflow  velocities  for  the  ro¬ 
tor  and  the  stator.  uZTr  and  u„,  are  the  axial  self-induced 
velocities  of  the  individual  components,  u and  uIjr  are 
the  axial  velocity  induced  by  the  stator  on  the  rotor  and 
v  ice  versa.  Xr  is  the  local  radius  of  the  individual  com¬ 
ponents.  X Rhr  and  Xrs,  are  the  radii  of  the  rotor  and 
stator  at  the  hub. 

To  determine  the  final  stator  diameter,  It  is  required 
that  the  tip  vortices  generated  by  the  rotor  blades  shall 
not  impinge  on  the  stator  blades.  In  other  words,  there 
should  not  be  large  velocity  gradients  over  the  tips  of  the 
stator  blades  so  that  the  stator  will  not  operate  in  the  tip 
vortices  from  the  rotor.  To  be  on  the  safe  side,  the  final 
stator  diameter  should  be  selected  so  as  to  be  smaller 
than  the  preliminary  stator  diameter  determined  based 
on  mass  conservation. 

Circulation  conservation  determines  the  magnitude  of 
the  stator  circulation  once  the  magnitude  of  the  rotor 
circulation  is  specified.  In  other  words,  the  magnitude  of 
the  stator  circulation  has  to  be  calculated  such  that  the 
total  circulation  is  conserved. 

DFSIGN  PROCEDURE 

A  flow  chart  of  the  design  method  for  PS  propulsors, 
which  is  similar  to  that  for  contrarotating  propellers  (see 
Chen  and  Reed  [4]),  is  shown  in  Figure  1.  It  consists  of 
three  phases:  specification  of  operating  conditions,  de- 
sien  and  analysis. 


(2) 

(3) 


Fig.  1.  Flow  chart  of  PS  propulsor  design  method. 


Specification  of  Operating  Conditions 

In  the  first  phase,  the  design  requirements  and  the 
wake  survey  data  need  to  be  provided.  The  effects  of 
the  hull  on  the  flow  and  hull-propulsor  interaction  are 
traditionally  represented  by  the  nominal  wake  and  two 
interaction  coefficients;  the  thrust  deduction  factor  de¬ 
scribed  previously  and  the  wake  fraction. 

The  wake  fraction  results  from  having  the  propulsor 
operate  in  the  spatially  non-uniform  wake  field  of  the 
hull.  This  non-uniform  wake  mainly  results  from  the 
shaft  inclination  and  boundary  layer  growth  along  the 
hull.  The  same  propulsor  model  can  be  operated  in  uni¬ 
form  flow  at  the  same  angular  velocity  and  the  advance 
speed  varied  until  the  thrust  agrees  with  that  measured 
behind  the  hull.  The  difference  between  the  ship  and 
thrust- identity  speeds  can  be  expressed  as  the  Taylor 
wake  fraction 


(4) 
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where  VA  is  the  advance  speed  and  V,  the  ship  speed. 


xrr 


The  wake  at  the  location  of  the  propulsor  in  the  ab¬ 
sence  of  the  propulsor  is  called  the  nominal  wake.  The 
effective  wake  distribution  determined  in  the  presence  of 
the  propulsor  depends  on  the  mutual  interaction  of  the 
propulsor  and  the  stern  boundary  layer.  In  the  design 
of  a  wake-adapted  propulsor,  the  radial  distribution  of 
effective  wake  is  traditionally  represented  by 

(1  -  wx(XR))  =  i^I(l  -  wc(XR))  (5) 

1  —  u>„ 

where  (1  —  u;c(A'r))  is  the  nondimensional  circumferential 
mean  axial  velocity  component  (nominal  wake)  from  the 
wake  survey,  1  —  wt  is  the  Taylor  wake  and  1  —  wv  the 
volumetric  mean  nominal  wake. 

Huang  et  al.[5]  developed  a  theoretical  calculation  of 
the  effective  wake  by  solving  a  simplified  Euler  equation 
in  a  conservative  force  field.  Moreover,  Shih  [6]  solved 
the  Euler  equation  in  a  non-conservative  force  field  to 
predict  the  velocity  distribution  of  a  propeller  operat¬ 
ing  in  an  inviscid,  axisymmetric  shear  flow.  Results  of 
both  theoretical  predictions  from  Huang  et  al.  and  Shih 
show  good  agreement  with  experimental  data.  However, 
the  calculation  of  the  effective  wake  for  PS  propulsors  is 
still  beyond  the  state-of-the-art  due  to  the  complex  in¬ 
teractions  between  the  PS  propulsor  and  the  hull,  and 
between  the  rotor  and  the  stator.  In  the  present  design 
procedure,  the  method  for  calculating  the  effective  wake 
for  a  SR  propeller  is  employed. 

Design 


u 


(b)  Stator. 


In  the  design  phase,  there  are  three  design  stages;  pre¬ 
liminary,  intermediate,  and  final. 

Preliminary  Design 

In  the  preliminary  design  stage,  one  investigates  a  lim¬ 
ited  number  of  design  parameters  such  as  diameter,  an¬ 
gular  velocity,  number  of  blades  and  radial  distribution  of 
loading.  In  general,  lifting-line  theory  can  economically 
determine  an  optimum  design  through  a  large  number 
of  parametric  studies  from  the  calculations  of  efficiency, 
strength,  and  cavitation.  Economy  of  effort  is  especially 
crucial  for  the  PS  propulsor  design  because  of  the  large 
number  parameters  as  well  as  the  calculation  of  mutual 
interactions  between  the  rotor  and  th<  stator. 

Kerwin  et  al.  [7]  developed  an  advanced  lifting-line 
theory  for  PS  propulsors.  A  vortex  lattice  method  was 
used  to  solve  the  lifting-line  theory.  The  induced  veloc¬ 
ities  in  the  theory  include  the  self-induced  velocities  of 
each  blade  row  as  well  as  the  velocities  induced  by  the 
other  blade  row.  A  velocity  diagram  for  PS  propulsors  is 
given  in  Figure  2. 


Fig.  2.  Velocity  diagram  for  PS  propulsors. 

Some  highlights  of  Kerwin  et  al.’s  theory  are  described 
in  the  following.  Three  fundamental  assumptions  are 
adopted: 

1.  The  blades  of  each  blade  row  axe  represented  by 
straight,  radial  lifting  lines. 

2.  The  blades  of  each  blade  row  are  considered  to  have 
equal  angular  spacing  and  identical  loading. 

3.  The  wake  geometry  is  assumed  to  be  purely  heli¬ 
cal,  with  the  pitch  at  each  radius  determined  either 
by  the  undisturbed  inflow  (linear  or  lightly  loaded 
theory),  or  by  the  induced  flow  at  the  lifting  line 
(moderately  loaded  theory). 

Based  on  a  thrust  identity  (it  can  also  be  represented 
by  a  torque  identity),  Kerwin  et  al.  obtained  the  opti¬ 
mum  circulation  using  a  variational  approach.  In  other 


words,  for  a  prescribed  total  thrust,  T,  and  a  prescribed 
torque  ratio,  q,  between  the  two  blade  rows,  the  opti¬ 
mum  circulation  of  each  blade  row  can  be  obtained  by 
minimizing  the  total  power. 

I  he  variational  approach  develops  as  follows.  Blade 
row  1  represents  the  rotor  and  blade  row  2  represents  the 

stator 

Minimize 

^iQi  +  W2Q2,  (6) 

subject  to  the  constraints: 

T,  +  T2  -  T  =  0.  (7) 

and 

<?2  -?<?!=  0.  (8) 

A  constrained  objective  function  is  formed: 


II  — !iQ\  —  — ' zQ :  +  A]  (Tj  +  7A  —  T)  +  A2( ^2  —  <lQ\ )i  (9) 


where  uri7r-  J„  and  u„m,jn  stand  for  the  influence  functions 
for  the  axial  and  the  tangential  velocities  induced  on  the 
control  point,  m,  of  the  ith  component  by  a  nth  element  of 
unit  strength  on  the  jth  component.  The  influence  func¬ 
tions  represent  the  self-induced  velocities  when  i  equals 
j  and  the  mutual  interaction  velocities  wnen  i  does  not 
equal  j. 

The  constrained  objective  function  H,  Eq  (9),  is  ex¬ 
panded  from  Eqs.  (10),  (11),  (12),  and  (13)  and  its  par¬ 
tial  derivatives  with  respect  to  the  unknown  circulations, 
r„,  and  Lagrange  multipliers,  A,,  are  set  equal  to  zero. 
This  process  provides  Mj  +  A/2  -I-  2  equations  which  can 
be  solved  for  the  circulations  at  the  control  points  and 
the  Lagrange  multipliers.  That  is  to  say.  one  obtains 

=  0  (1  =  1 . A/,,j  =  1,2).  (14) 

and 

^=0  (1  =  1,2).  (15) 


where  —q  and  -2  are  the  angular  velocities  for  each  blade 
row,  Aj  and  .\3  are  Lagrange  multipliers,  7)  and  T3  axe  the 
thrusts  for  each  blade  row,  and  Qx  and  are  the  torques 
tor  each  blade  row.  Thrusts  and  torques  are  computed 
using  the  following  formulas: 

M, 

~  p  y  ■  ZA  im  -A  A  (  V(,m  -f-W,  X  ntm  Tu<im  )  (:  =  1,2), 

n;  =  1 

(10) 

and 

M, 

f-1  /  .  rirn  A  A  nxm  A  f{XTn  ( \  rjm  -f-  urim )  (z  =  1 , 2). 

m  —  i 

(id 

where  X n,m  is  the  local  radius  from  the  shaft  axis  to  the 
control  point,  m,  XXn,m  is  the  radial  distance  between 
the  two  lattice  points  surrounding  control  points  m,  Z, 
is  the  number  of  blades  and  M,  the  number  of  control 
points  r,„,  is  the  circulation,  Vxlm  and  V,,m  are  the  axial 
arid  the  tangei  il  inflow  velocities,  and  uxlm  and  u(lm 
are  the  axial  aim  the  tangential  induced  velocities  at  the 
control  points.  These  induced  velocities  include  the  self- 
mducc'd  velocities  and  the  velocities  induced  by  the  other 
blade  row.  The  velocities  uxim  and  u„m  can  be  expressed 
as  follows. 

-  Ah 

Ur.™  =  rjnuximjn  (f  =  1,2),  (i2) 


7  , 

(z  =  1,2),  (13) 

; = 1 n= 1 


The  optimum  circulation  distribution  in  Kerwin  et 
al.’s  theory  allows  non- zero  loading  at  the  blade  roots 
which  is  an  important  factor  for  PS  propulsors.  The  the¬ 
ory  provides  an  equal  and  opposite  circulation  for  the 
rotor  and  the  stator  at  the  hub  to  ensure  minimum  hub 
vortex  strength.  In  other  words,  the  net  circulation  for 
the  two  individual  blade  rows  at  the  hub  is  zero.  There¬ 
fore,  the  pr  jsent  design  procedure  accounts  for  the  effects 
of  the  hub  boundaries.  In  addition,  the  slope  of  the  cir¬ 
culation,  is  constrained  to  be  essentially  zero  at  the 
hub  so  that  the  trailing  vortex  sheet  in  the  hub  region 
may  be  ’’eliminated”. 

The  optimum  circulations  for  a  SR  propeller  between 
the  Kerwin  et  al.  and  the  Caster  et  al.[Sj's  programs  were 
compared  in  Chen  et  al.[9].  The  results  from  Kerwin 
et  al.’s  program  are  very  good.  This  indicates  that  the 
variational  method  is  an  appropriate  approach  to  solve 
lifting-line  theory.  However,  this  does  not  indicate  that 
Kerwin  et  al.’s  program  is  an  adequate  tool  to  design  PS 
propulsors  unless  it  can  be  validated.  This  validation  is 
part  of  the  effort  pursued  in  the  present  study. 

The  real  flow  passing  through  the  rotor  and  the  stator 
can  be  envisioned  as  an  unsteady  flow  both  in  time  and 
space  because  of  the  different  blade  numbers  and  rotation 
speeds.  However,  as  is  usually  done,  the  design  problem 
is  simplified  by  the  use  of  circumferential-average  inflow. 
Calculating  the  mutual  interaction  velocities  is  one  of  the 
most  critical  factors  affecting  the  PS  propulsor  design.  To 
obtain  accurate  mutual  interaction  velocities,  the  use  of 
the  interaction  velocities  calculated  from  lifting-surface 
theory  developed  by  Chen  and  Reed  [10]  is  employed  in 
the  present  study. 
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Intermediate  Design 

In  the  intermediate  design  stage,  cavitation  and  strength 
are  major  concerns.  Blade  surface  cavitation  results  in 
blade  erosion,  noise,  and  thrust  loss  which  are  detrimen¬ 
tal  to  ship  performance.  In  order  to  improve  the  cavi¬ 
tation  performance,  one  can  vary  blade  thickness,  chord 
distribution,  and  blade  loading.  However,  the  strength 
and  the  propulsive  performance  are  also  affected  by  vary¬ 
ing  the  above  parameters. 

As  the  propeller  rotates,  it  is  subjected  to  both  hy¬ 
drodynamic  and  centrifugal  loadings.  To  keep  the  stress 
within  a  blade  below  a  certain  allowable  level,  a  propeller 
blade  must  contain  enough  material.  This  acceDtable 
stress  level  is  affected  by  the  material  properties  which 
are  functions  of  steady  state  loading,  fatigue  strength, 
mean  and  unsteady  blade  loadings.  The  propeller  geo¬ 
metric  parameters,  such  as  chord  length,  thickness,  skew, 
and  rake  will  affect  the  stress.  The  stress  can  be  calcu¬ 
lated  approximately  using  simple  beam  theory  which  rep¬ 
resents  the  propeller  blade  as  a  straight  cantilever  beam 
with  variable  cross  section  without  camber.  However,  a 
finite  element  analysis  is  required  for  a  more  accurate 
analysis  for  the  full  power  ahead  and  the  backing  cases. 

Even  for  SR  propellers,  it  is  well  recognized  that  cav¬ 
itation  predictions  frequently  do  not  match  the  measure¬ 
ment  well.  This  discrepancy  is  mainly  due  to  the  lim¬ 
itations  of  the  prediction  tools.  Cavitation  prediction 
is  especially  complicated  in  the  PS  propulsor  design  be¬ 
cause  of  the  interactions  between  the  rotor  and  the  stator. 
Owing  to  the  temporal  and  spatial  variations  in  the  ve¬ 
locity  field  of  the  rotor  and  the  stator,  a  method  which 
can  handle  the  circumferential  variations  in  the  mutual 
interaction  velocities  for  PS  propulsors  has  to  be  devel¬ 
oped.  This  is  beyond  the  current  state-of-art. 

Final  Design 

In  the  final  design  stage,  a  lifting-surface  theory  which 
incorporates  three  -dimensional  flow  field  effects  was  em¬ 
ployed  to  determine  the  detailed  blade  geometry  (pitch 
and  camber  distributions).  In  the  present  study,  a  lifting- 
surface  program  with  hub  image  effects  developed  by 
Wang  [11]  is  employed.  Blades  and  their  wakes  are  rep¬ 
resented  by  vortex  lattices.  The  hub  is  represented  by 
a  distribution  of  dipoles  which  ends  at  the  hub  apex.  A 
non-zero  circulation  at  the  hub  of  a  propulsor  is  set  so 
that  the  normal  velocity  at  the  hub  will  be  zero. 

Though  the  effect  of  the  viscous  drag  has  been  consid¬ 
ered  in  the  force  calculations,  all  the  design  codes  used 
in  this  effort  are  based  on  potential  flow  theory.  The 
viscous  drag  has  been  calculated  using  the  following  em¬ 
pirical  formula. 

CD  =  C,{  1  +  1.25  (t/c)  +  125(t/c)“)  (16) 


where  C/  is  the  skin  friction  coefficient  for  a  smooth  plate. 


and  t/c  is  the  section  thickness  to  chord  ratio.  The  value 
of  Cj  varies  from  0.004  to  0.008  for  corresponding  blade 
section  Reynolds  number  varying  between  108  and  106. 

In  the  design  procedure,  boundary  layer  calculations 
need  to  be  performed  to  ensure  unseparated  flow  along 
the  blade.  A  program  developed  by  Cebeci  [12]  was  em¬ 
ployed  in  the  present  study.  This  program  can  calcu¬ 
late  incompressible  laminar  and  turbulent  boundary  lay¬ 
ers  on  plane  and  axisymmetric  bodies  with  either  smooth 
or  rough  surfaces. 

Analysis 

In  the  analysis  phase,  steady  forces  and  unsteady  forces 
and  moments  need  to  be  calculated  using  inverse  lifting- 
surface  codes.  To  determine  the  resultant  steady  thrust, 
torque,  and  efficiency  of  the  PS  propulsor  under  design 
and  off-design  conditions,  the  vortex  lattice  method  in¬ 
cluding  hub  effects,  developed  by  Greeley  and  Kerwin 
[13]  was  employed. 

Fluctuating  forces  and  moments  arise  from  wake  non¬ 
uniformities  or  shaft  inclination.  These  unsteady  compo¬ 
nents  can  cause  cavitation,  strength,  and  radiated  noise 
problems.  In  general,  larger  skew  produces  smaller  vibra¬ 
tory  forces  and  moments,  while  producing  larger  blade 
stresses.  Hence,  from  considerations  of  both  vihratorv 
forces  and  blade  stresses,  the  skew  distribution  needs  to 
be  optimized.  The  design  is  complete  when  the  unsteady 
shaft  forces  and  moments  are  below  the  design  require¬ 
ments.  A  method  which  can  handle  the  variations  in  the 
mutual  interaction  velocities  between  the  rotor  and  the 
stator  is  the  key  point  in  calculating  the  unsteady  forces 
and  moments  for  PS  propulsors. 


A  POSTSWIRL  PROPULSOR  DESIGN 

A  PS  propulsor  design  is  presented  using  the  method 
developed  in  the  previous  section.  Some  of  the  design 
stages  will  be  skipped  because  the  present  design  is  for 
close- to-uniform  flow. 

DESIGN  REQUIREMENTS 

The  PS  propulsor  was  designed  for  close-to-umiorm 
flow  at  the  operating  point  for  a  surface  ship.  The  ship 
speed  was  chosen  as  20  knots  (10.3  m/s).  The  thrust 
loading  coefficient,  Cjk,  is  0.2871.  The  rotor  diameter  is 
7.55  ft  (2.30  m)  and  the  rotation  speed  199  rpm.  The 
blade  numbers  of  the  rotor  and  stator  are  5  and  7,  re¬ 
spectively.  The  ship  full  power  condition  is  2,970  shaft 
horsepower  (2,216  KW)  per  shaft. 

The  design  parameters  for  the  present  study  were  cho¬ 
sen  based  on  a  parametric  study.  The  stator  diameter 
was  determined  through  mass  conservation,  in  Equation 
(2).  To  ensure  that  the  stator  operates  inside  the  tip  vor- 
ticr  of  the  rotor,  the  final  stator  diameter,  which  is  85 
percent  of  the  rotor  diameter,  was  chosen  to  be  slightly 
smaller  than  the  preliminary  diameter  calculated  using 


5 


mass  conservation.  The  axial  spacing  was  chosen  to  be 
one  quarter  of  the  rotor  diameter.  A  summary  of  the 
design  parameters  for  the  PS  propulsor  is  given  in  Table 
1. 


Table  1.  PS  propulsor  design  —  Summary. 


Rotor 

Stator 

Ship  Speed  (knots) 

20 

56 

Rotational  Speed  (rpm) 

199 

0 

Thrust  Loading 

0.2871 

Coefficient 

Geometry 

Diameter  (ft) 

7.55 

6.42 

Number  of  blades 

[  5 

7 

Expanded  Area  Ratio 

0.538 

0.611 

Skew  (deg) 

25 

0 

Total  Rake 

0 

0 

Blade  Sections 

* 

* 

Axial  Spacing  (ft) 

\ _ _ 

*  NACA  66  (TMB  Modified)  Thickness, 
NACA  a  =  0.8  Meanline 


PRELIMINARY  DESIGN 

The  lifting-line  calculations  will  be  discussed  in  this 
section.  The  optimum  circulation  distributions  for  the 
rotor  and  the  stator  are  determined  by  the  variational  ap¬ 
proach  described  in  the  previous  section.  The  optimum 
and  the  root-  and  tip-unloaded  circulation  distributions 
for  the  rotor  and  stator  are  shown  in  Figure  3.  The  ad¬ 
vantages  of  unloading  the  blade  root  and  tip  are  to  delay 
the  blade  root  and  tip  cavitation  inception  and  to  reduce 
the  tendency  toward  cavitation  erosion  near  the  blade 
root  and  tip.  Since  the  lift  coefficient  at  the  blade  root 
was  not  very  high  with  the  optimum  circulation  distribu¬ 
tion,  a  20  percent  reduction  in  the  blade  loading  at  the 
root  was  made.  In  addition,  the  slope  of  the  circulation, 
was  constrained  to  be  essentially  zero  at  the  blade 
root  to  minimize  the  trailing  vortex  sheet. 

The  guidelines  for  unloading  the  stator  hub  have  to 
consider  two  factors:  (1)  ensure  zero  hub  vortex  strength 
by  providing  an  equal  and  opposite  total  circulation  at 
the  hub  with  respect  to  the  rotor,  and  (2)  minimize  the 
trailing  vortex  sheet  by  maintaining  almost  zero  circula¬ 
tion  slope  at  the  hub. 

Based  on  cavitation,  flow  separation,  and  efficiency 
considerations,  the  chord-length  distributions  of  the  rotor 
and  stator  were  chosen.  The  thickness  distribution  was 
selected  based  on  strength  and  cavitation  considerations. 

A  skew  distribution  of  25  degrees  at  the  tip,  varying 
nonlinearly  from  zero  at  the  hub,  was  selected  for  the 
rotor.  The  stator  was  designed  with  zero  skew.  The 
total  rake  for  both  the  rotor  and  the  stator  was  zero. 
Therefore,  the  rotor  was  given  negative  rake  to  offset  the 
skew-induced  rake. 


Nondimensional  radius,  Xr 


Fig.  3.  Optimum  and  unloaded  circulation 
distribution  for  rotor  and  stator. 


INTERMEDIATE  DESIGN 

In  the  intermediate  design  stage,  strength  and  cavi¬ 
tation  was  considered.  The  stresc  distributions  for  the 
rotor  and  stator  were  calculated  and  were  determined  to 
be  well  below  an  allowable  stress  of  12,500  psi  maximum 
stress  for  Nickel-Aluminum-Bronze  material.  Although 
this  propulsor  set  was  designed  for  close-to- uniform  flow, 
the  cavitation  performance  of  the  stator  is  important 
due  to  the  non-uniform  effects  of  the  rotor  on  the  sta¬ 
tor.  Since  there  is  no  appropriate  method  to  predict  the 
cavitation  for  the  stator,  the  cavitation  prediction  was 
not  carried  out  for  this  design. 

FINAL  DESIGN 

The  final  pitch  and  camber  distributions  were  deter¬ 
mined  using  lifting-surface  theory  with  hub  effects  in¬ 
cluded.  For  marine  propeller  applications,  an  a  =  0.8 
meanline  loading  is  an  appropriate  chordwise  loading  dis¬ 
tribution  from  cavitation  and  viscous  flow  points  of  view. 
A  NACA  66  thickness  distribution  was  selected  for  the 
present  application. 

The  final  rotor  camber  and  pitch  distributions  are 
shown  in  Figures  4  and  5  while  the  final  stator  camber 
and  hydrodynamic  pitch  angle  distributions  are  shown  in 
Figures  4  and  6. 


Maximum  camber/Chord  length,  fvi/c 


Nondimensional  radius,  Xr 


Fig.  4.  Final  maximum  camber  to  chord  length 
ratio  distributions. 


Nondimensional  radius.  Xr 


Fig.  5.  Final  pitch  to  diameter  ratio  for  rotor. 


Nondimensional  radius,  Xr 


Fig.  6.  final  hydrodynamic  pitch  angle 
for  stator. 


PERFORMANCE  PREDICTION  AND 
EXPERIMENTAL  RESULTS 

Aluminum  models  of  the  rotor  and  stator  were  man¬ 
ufactured  based  on  the  final  design  geometry.  The  sell- 
propulsion  tests  were  carried  out  using  DTRC  model  pro¬ 
peller  5118  to  represent  the  rotor  and  model  propeller 
5119  to  represent  the  stator.  The  self-propulsion  tests 
were  performed  at  DTRC. 

SELF-PROPULSION  TESTS 

The  resistance  and  the  self-propulsion  tests  for  the  de¬ 
signed  PS  propulsor  were  performed  using  the  same  ship 
Model  5365-A  as  those  for  the  stock  SR  propeller.  Figure 
7  shows  the  measured  delivered  power  and  rotation  speed 
of  the  unit  as  a  function  of  ship  speed  as  well  as  the  pre¬ 
dicted  performance  at  the  design  ship  speed.  Table  2  also 
shows  the  predicted  and  measured  self-propulsion  perfor¬ 
mance  at  the  design  ship  speed  of  the  PS  propulsor. 

The  effective  horsepower  used  in  the  design  is  4.6  per¬ 
cent  lower  than  the  measurement.  This  difference,  due 
to  the  modification  of  the  test  vehicle,  may  cause  some 
discrepancy  between  the  design  and  the  experiment.  The 
delivered  power  from  the  prediction  is  1 .3  percent  higher 
than  that  from  the  experimental  values.  The  thrust  de¬ 
duction  and  the  wake  fraction  estimated  for  the  design 
are  very  close  to  the  experimental  values. 
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Rotor  rotation  speed,  (rpm)  Delivered  power,  PD  (hp) 


Ship  speed,  Vs  (knots) 
(a)  Delivered  power. 


Ship  speed,  V$  (knots) 

(b)  Rotation  speed. 

Fig.  7.  Predicted  and  measured  self-propulsion 
test  results. 


Table  2.  Predicted  and  measured  powering 

performance  of  PS  propulsor  design. 


Prediction 

Self-propulsion  Experiment 

Vs  (knots) 

20 

20 

PH(hp) 

1,550  (-4.6%) 

1,625 

Pd  (hp) 

2,299  (  1.3%) 

2,269 

1-t 

0.865  (-0.6%) 

0.870 

I-wx 

1.000  (-0.5%) 

1.005 

Ja 

1.349  (-1.5%) 

1.370 

N  (rpm) 

199  (  1.0%) 

197 

Kt 

0.207  (-5.0%) 

0.218 

Kq 

0.056  (-1.8%) 

0.057 

T)D 

0.674  (-5.7%) 

0.715 

0.788  (-4.5%) 

0.825 

TlR 

1.000 

1.000 

The  design  advance  coefficient  is  1.5  percent  lower  and 
the  predicted  rotation  speed  is  1.0  percent  higher  than 
the  measurement.  The  thrust  and  the  torque  predictions 
are  5.0  percent  and  1.8  percent  lower  than  the  experi¬ 
ment.  The  predicted  open  water  efficiency  and  propul¬ 
sive  efficiency  are  4.5  percent  and  5.7  percent  lower  than 
the  measurement.  These  discrepancies  are  primarily  due 
to  the  difference  between  the  design  and  the  measured 
effective  horsepower.  In  general,  the  accuracy  of  the  ex¬ 
perimental  measurements  with  the  PS  propulsor  is  ±  2 
percent  on  thrust  and  torque.  Overall  the  predicted  val¬ 
ues  agree  well  with  the  experimental  measurements,  and 
in  general  are  within  the  accepted  accuracy  of  the  ex¬ 
perimental  measurements.  A  view  of  the  PS  propulsor 
created  on  a  Computer  Vision  system  is  shown  in  Figure 
8. 
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Fig.  8.  A  configuration  of  PS  propulsor. 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  following  conclusions  can  be  drawn  from  the  present 

study. 

1  The  present  design  method  is  an  adequate  proce¬ 
dure  for  the  PS  propulsor  design.  Self-propulsion 
experiments  show  that  the  performance  predictions 
agree  well  with  the  experimental  measurements. 

2.  The  mutual  interaction  velocities  between  the  ro¬ 
tor  and  the  stator  play  important  roles  in  the  PS 
propulsor  design.  This  is  particularly  important  for 

There  are  two  general  recommendations  for  future  work. 

1.  Propulsor  -  hull  interactions  and  effective  wake:  A 
method  for  calculating  the  PS  propulsor  -  hull  inter¬ 
actions,  including  thrust  deduction  and  wake  frac¬ 
tion,  and  effective  wake  needs  to  be  developed.  This 
is  particularly  important  for  the  stator  because  the 
complex  interactions  between  the  PS  propulsor  and 
the  hull,  and  the  rotor  and  the  stator. 

2.  Cavitation  and  unsteady  force:  A  method  which  can 
handle  the  circumferential  variations  of  the  mutual 
interaction  velocities  between  the  rotor  and  the  sta¬ 
tor  needs  to  be  developed.  This  method,  as  a  mat¬ 
ter  of  fact,  is  crucial  to  the  cavitation  and  unsteady 
force  calculations  for  the  stator. 
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APPENDIX 


Parametric  Studies 

The  purpose  of  this  study  is  to  secure  an  optimum  design  for  a  PS  propulsor  through 
parametric  studies.  The  lifting-line  theory  developed  by  Kerwin  et  al.  as  in  reference  (7) 
was  employed.  The  parametric  studies  covered  the  following  range  of  parameters. 

•  Rotor  diameter,  Dr ,  6.5  -  8.0  ft  (1.98  -  2.44  m). 

•  Rotor  rotation  speed,  Nr,  150  -  210  rpm. 

•  Ratio  of  stator  and  rotor  diameter,  Da/Dr,  0.85  -  0.95 

As  shown  in  Figure  9,  although  the  selected  rotor  diameter,  7.55  ft  (2.30  m),  is  not  the 
optimum  diameter,  the  difference  of  the  propulsive  efficiency  between  the  optimum  and  the 
selected  diameters  is  around  0.3  point. 

Figure  10  shows  that  the  selected  rotor  rotation  speed,  199  rpm,  is  not  the  optimum 
rotation  speed.  However,  the  difference  of  the  propulsive  efficiency  between  the  optimum 
and  the  selected  rotation  speeds  is  small. 

The  optimum  rotor  diameter  and  rotor  rotation  speed  were  calculated  only  based  on  the 
consideration  of  efficiency.  However,  the  selected  rotor  and  rotation  speed  were  calculated 
based  on  consideration  not  only  of  efficiency  but  of  cavitation,  strength,  and  flow  separation. 
Therefore,  there  is  a  slight  difference  between  the  optimum  and  selected  rotor  diameter  and 
rotation  speed. 

Figure  11  indicates  that  the  propulsive  efficiency  is  insensitive  to  the  ratio  of  the  sta¬ 
tor  and  rotor  diameters.  Since  lifting-line  theory  does  not  account  for  mass  conservation 
between  the  rotor  and  stator,  additional  consideration  must  be  included. 

The  basic  guideline  for  determining  the  final  stator  diameter  is  to  ensure  that  the  rotor 
tip  vortices  will  not  impinge  on  the  stator  blades  because  they  can  result  in  erosion  or  blade 
rate  noise  of  the  stator.  To  be  on  the  safe  side,  the  final  stator  diameter,  was  selected  to  be 
85  percent  of  the  rotor  diameter.  This  is  slightly  smaller  than  the  rotor  slipstream  diameter 
calculated  using  mass  conservation.  As  seen  in  Figure  11,  there  is  almost  no  change  in 
efficiency  with  varying  the  ratio  of  the  stator  diameter  to  the  rotor  diameter.  In  other 
words,  there  is  no  penalty  in  efficiency  to  select  this  specific  stator  diameter.  The  axial 
spacing  was  chosen  to  be  one  quarter  of  the  rotor  diameter. 


Geometric  Specification 

Based  on  cavitation,  flow  separation,  and  efficiency  considerations,  the  chord-length 
distributions  of  the  rotor  and  stator  were  chosen.  The  chord  distributions  for  the  rotor  and 
stator  are  shown  in  Figure  12.  The  expanded  area  ratio  (EAR)  was  calculated  based  on 
both  Burrill’s  and  Keller’s  criteria  for  cavitation  inception. 

The  thickness  distribution  was  selected  based  on  strength  and  cavitation  considerations. 
Figure  13  shows  the  thickness  distributions  for  the  rotor  and  stator. 

As  shown  in  Figure  14  ,  a  tip  skew  distribution  of  25  degrees,  varying  nonlinearly  from 
zero  at  the  hub,  was  selected  for  the  rotor.  The  stator  was  designed  with  zero  skew.  The 
total  rake  for  both  rotor  and  stator  was  zero.  Therefore,  the  rotor  has  negative  rake  to  offset 
the  skew-induced  rake.  The  stress  distributions  computed  by  beam  theory  corresponding 
to  these  choices  of  geometry  are  given  in  Figure  15. 
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Fig.  11.  Propulsive  efficiency  versus  ratio  of  stator 
and  rotor  diameter  for  fixed  rotor  diameter 
and  rotation  speed. 
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Fig.  12.  Radial  distribution  of  chord  length. 
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Fig.  13.  Radial  distribution  of  thickness. 
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Fig.  14.  Radial  distribution  of  skew. 
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Fig.  15.  Radial  distribution  of  stress. 
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Final  Postswirl  Propulsor  Geometry 

The  final  geometric  specifications  of  the  PS  propulsor,  including  the  details  of  the 
leading  and  ^ailing  edges  were  computed  using  the  computer  code,  XYZ-PROP,  developed 
by  Brackett.*  All  the  input  data,  such  as  the  chord  length,  thickness,  skew,  pitch  and 
camber  distributions,  were  faired  by  a  cubic  spline  procedure  before  being  input  to  XYZ- 
PROP.  A  list  of  the  chord  length,  thickness,  skew,  pitch  and  camber  distributions  are  given 
in  Table  3. 


Table  3.  Final  design  geometry  for  PS  propulsor. 

a.  Rotor. 


l  m — 

q/b 

P/D 

hvD 

_ 6s _ 

t/c 

WKusam 

TTfTTI 

■n 

■■CKiTiTtTi] 

■witioni 

0.22702 

0.01876 

0.2500 

0.16030 

0.0000 

0.67623 

0.19527 

0.02084 

0.3000 

0.17084 

1  tia 

0.0000 

1.18209 

0.17324 

0.02336 

0.4000 

0.19178 

0.0000 

2.31700 

0.13278 

0.02978 

0.5000 

0.21059 

0.0000 

3.76256 

0.09916 

0.03360 

0.6000 

0.22638 

1.7875 

0.0000 

5.87000 

0.07453 

0.03250 

0.7000 

0.23900 

1.6875 

0.0000 

8.99300 

0.05885 

0.02519 

0.8000 

0.24709 

0.0000 

13.31010 

0.05067 

0.01632 

0.9000 

0.23416 

0.0000 

18.82284 

0.04747 

0.00879 

0.9500 

0.19326 

0.0000 

21.86989 

0.04681 

0.00595 

1.0000 

0.00000 

1.1000 

0.0000 

25.00000 

0.04635 

0.00342 

b.  Stator. 


r/R 

c/D 

P/D 

iT/D 

0s 

t/c 

0.3281 

0.25156 

■a 

0.0000 

0.07459 

-0.01055 

^■ik  .in 

0.24893 

96.4790 

0.0000 

0.0000 

0.07328 

-0.01100 

0.24278 

95.7450 

0.0000 

0.0000 

0.07046 

-0.01478 

0.23647 

94.9375 

0.0000 

0.0000 

0.06811 

-0.02145 

0.23009 

94.3715 

0.0000 

0.0000 

0.06625 

-0.02674 

0.21751 

93.8499 

0.0000 

0.0000 

0.06380 

-0.02990 

0.7000 

0.20585 

93.6830 

0.0000 

0.0000 

0.06273 

-0.02596 

0.8000 

0.19524 

93.6510 

0.0000 

0.0000 

0.06249 

-0.01940 

0.9000 

0.17476 

94.0000 

0.0000 

0.0000 

0.06286 

-0.01053 

0.9500 

0.14502 

94.6250 

0.0000 

0.0000 

0.06324 

-0.00469 

1.0000 

0.03795 

95.9800 

0.0000 

0.0000 

0.06369 

0.00223 

*  Brackett,  T.  E.  1976.  Analytical  Specification  of  Propeller  Blade-Surface  Geometry. 
DTRC  Ship  Hydrodynamics  Department  Report,  DTRC/SHD-699-01. 
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